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On top of the mean-field analysis of a Bose-Einstein condensate, one typically applies the Bogoli¬ 
ubov theory to analyze quantum fluctuations of the excited modes. Therefore, one has to diagonalize 
the Bogoliubov Hamiltonian in a symplectic manner. In our article we investigate the topology of 
these Bogoliubov excitations in inversion-invariant systems of interacting bosons. We analyze how 
the condensate influences the topology of the Bogoliubov excitations. Analogously to the fermionic 
case, here we establish a symplectic extension of the polarization characterizing the topology of 
the Bogoliubov excitations and link it to the eigenvalues of the inversion operator at the inversion- 
invariant momenta. We also demonstrate an instructive but experimentally feasible example that 
this quantity is also related to edge states in the excitation spectrum. 
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I. INTRODUCTION 


Since the discovery of Bloch bands with nontrivial 
topological structure, the field of topological insulators 
and superconductors has been rapidly growing (TJ [2] . The 
most prominent example is the integer quantum Hall ef¬ 
fect, where one can link the Hall conductance of the 
ground state with the Chern number of the occupied 
bands [3 J . This strict quantization is due to the fermionic 
character of the particles, forcing all states within a band 
to be equally occupied. For this reason, a system consist¬ 
ing of bosons does not exhibit such quantized observables. 

In addition, there are topologically protected edge 
states as a consequence of the bulk-boundary relation [I] 
2 . Thus, for noninteracting particles the band structure 
determines the existence of edge states. Although nonin¬ 
teracting bosons condense in the lowest-energy mode, the 
wave function of the excited modes within a band can ex¬ 
hibit a topological structure ED IS]- Therefore, there can 
be edge modes in the excitation spectrum of noninteract¬ 
ing bosons. 

In recent years, there has been a great effort in creating 
nontrivial topological structures of fractional quantum 
Hall states in strongly interacting bosonic systems EHT3]. 
In contrast, the investigation of the topology of bosons 
in the weakly interacting regime has received insignif¬ 
icant attention. An interesting approach is discussed 
in Ref. Q0], where the edge states of the Su-Schrieffer- 
Heeger model become dynamically unstable by properly 
preparing the condensate in an excited transverse mode. 
In such a setup one considers the excitations above the 
Bose-Einstein condensate (BEC), which are effectively 
single-particle-like due to an expansion of the Hamilto¬ 
nian in orders of the condensate density [15]. The re¬ 
sulting Bogoliubov Hamiltonian couples particle excita¬ 
tions and hole excitations and has to be diagonalized in 
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a symplectic manner due to the bosonic commutation 
relation. Therefore, the definition of a topological in¬ 
variant for these Bogoliubov-Bloch bands is a priori not 
clear. As a result, the condensed part of the atoms has a 
substantial influence on the topology of the Bogoliubov 
excitations, which has not been discussed in the previous 
literature. There are only few articles about the defini¬ 
tion of a Chern number for bosonic Bogoliubov bands, 
though in the context of magnonic systems 

In contrast, here we focus on the treatment of the 
topology in inversion-invariant systems of weakly inter¬ 
acting bosons in one dimension. Our findings are relevant 
for condensed matter simulations with cold atoms in op¬ 
tical lattices fl8ff22] . 

In fermionic systems, the topological invariant of 
inversion-invariant systems is given by the macroscopic 
polarization [23] [24], which is a geometric phase of the oc¬ 
cupied bands [25]. Although its strict quantization is due 
to the fact that all orbitals of the bands below the Fermi 
energy are equally occupied, one can also consider it to 
be a quantity describing the structure of the Bloch bands 
independent of the fermionic character of the particles. 
For this reason, for bosonic systems it can be considered 
to be a topological invariant not directly connected to 
a bulk observable, but predicting the existence of edge 
states. 


Here we treat the bosonic condensate within a mean- 
field approach. In this context we note that Refs. [26H29] 
study the topology of fermionic systems using mean-field 
approximations. 


The article is organized as follows. In Sec. |TT| we in¬ 
troduce a model that is considered as an instructive but 
experimentally feasible example throughout the article. 
This system has the property that the condensate in¬ 
fluences the topology of its Bogoliubov excitations. In 


Sec. Ill we recall Bogoliubov theory. In Sec. IV we in¬ 
vestigate the topology of the Bogoliubov excitations. A 
main result of our article is Eq. (20), which defines an ex¬ 
tension of the macroscopic polarization for bosonic Bo¬ 
goliubov excitations. We also discuss the problems for 
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FIG. 1: (Color online) (a) Sketch of the Hamiltonian 0. 
(b) Corresponding phase diagram distinguishing between dif¬ 
ferent wave functions of the condensate for the parameters 
Vso/w — 0.2 and Xs,s' — ^s,a'X- ( c ) Instances of the corre¬ 
sponding wave functions for vq/uj = 0.1. The labels I, II, and 
III denote phases with a localized wave function, condensation 
at zero momentum, and condensation at finite momentum 
k > 0, respectively, (d) For phase II, we depict some Bogoli- 
ubov dispersion relations. The corresponding parameters are 
marked with arrows in the phase diagram. As explained in the 
main text, one can distinguish the topology of the Bogoliubov 
excitations depending on the symplectic polarization P s . For 
this reason, one can split phase II into a trivial phase Ha and 
a topological phase lib. The phase boundary is strongly influ¬ 
enced by the condensate due to the interactions between the 
particles. In Sec. GYS we show that topologically protected 
edge states appear in the topological phase lib. 


In position space the Hamiltonian reads 

M/2 

H= H m + V m , (1) 

m = —M/2+1 


where 


Hjn - cu 

— uq (<^+ m a+ m+ i — n| ?m n^ 5m+ i + H.c.^ 

^so 1 l,rn+l H.C.J , 


Vm = 


s,s / =‘t,t 




,m- 


( 2 ) 

( 3 ) 


Here m denotes the position of the wells and t, i denote 
the internal degree of freedom of the bosons. The two 
states have a level splitting of uj. The atoms can jump 
between the wells, described by the parameters uq and 
v so . The modulus of the hopping integral z/q is equal 
for the two spin components but differs in sign. The 
term proportional to v so denotes a spin-orbit coupling, 
which can be generated within current technology m 
ng. Additionally, we have a state-dependent interaction 
that is local in position space. Therefore, we ensure that 
this is not in conflict with the inversion symmetry of the 
Hamiltonian. 

The single-particle Hamiltonian corresponds to the 
systems in [311132] discussing fermionic systems of cold 
atoms that have a nontrivial topology. These articles 
suggest possible experimental implementations for the 
single-particle Hamiltonian. These could be also applied 
to bosonic systems. Additionally, one can control the 
interactions between the particles using Feshbach reso¬ 
nances. 


B. Mean-field expansion 

We shift the operators 


the definition of the polarization caused by the Gold- 


stone mode appearing in the lowest band. In Sec. IV E 


we apply our findings to the instructive model and show 
the existence of edge states. 


II. MODEL SYSTEM 


A,m 



Ct/4-,m + «f/4., 


,m •> 


( 4 ) 


where Ct/t,m €= C and Nq denotes the number of particles 
in the condensate which is assumed to be macroscopically 
occupied. The bosonic operators now account for quan¬ 
tum fluctuations on top of the condensate. We expand 
the Hamiltonian as 


A. Hamiltonian 

We consider systems of weakly interacting bosons in 
a periodic potential. An instance of such a system is 
sketched in Fig. [lja). There, we consider an ensemble 
of bosonic atoms with internal degree of freedom (spin) 
confined in an array of wells created by an optical lattice. 


H = pqEqp + y P^H W + H^ + O(po 1/2 ), (5) 

where /+> = Nq/M is the density of the condensed parti¬ 
cles. The Hamiltonian H^ (. H depends on {Ct/+m} 
and contains terms that are linear (quadratic) in bosonic 
operators. 

Here Eqp is a function of {Ct/f,m} an d denotes the 
Gross-Pitaevskii functional. It exactly reads as 0 with 
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the operators replaced by Ct/i,m an <4 Xs,s f Xs,s'Po- 
To find the mean-field ground state, we minimize it by 
an appropriate choice of The minimization proce¬ 

dure is performed by using a modified ansatz of Ref. [33] . 
Details can be found in Appendix |A| 

The result is depicted in the phase diagram in Fig.JTJb) 
for the special choice xtt = Xll = X and Xtt = 0- Note 
that for the mean-field treatment in Fig. [l we assume 
that nearly all bosons are condensed, so we approximate 
po ~ P = N/M , which is the density of all particles. 

We find three phases. In phase I appearing for poX < 0 
we find that the atoms condense within a small area 
with a localized wave function (34[ |35|. In phase II 
we find a condensation at zero momentum k = 0 and 
(Ct,m,Ct,m) = (1,0). Due to the spin-orbit coupling and 
the interactions, the atoms condense at a finite momen¬ 
tum k > 0 in phase III. Thus, only in phase II there 
is a mean-field wave function, which is invariant under 
inversion. 

At a stationary point of the Gross-Pitaevskii func¬ 
tional, the linear part in © vanishes and the excitations 
are solely governed by the quadratic Bogoliubov Hamil¬ 
tonian. To respect that we work at constant particle 
number we consider Nq to be an operator and replace m 

Nq = N — a Jn,s^rn,s- ( 6 ) 

m, s=td 


where a£ = is a vector of bosonic creation 

operators and 



— uj — cos k 2w so sin k 

—2iu so smk uj + 2 z^o cos k 


H« = 
H (2) = 


4 XP - Meff 0 

0 2 XUP - Meff 

2 XP 0^ 

0 0 ) • 


(9) 

( 10 ) 

( 11 ) 


The chemical potential reduces to /x e ff = — uo — 2 uq + 2py. 
The Bogoliubov Hamiltonian determines the topologi¬ 
cal properties of the excitations. In order to investigate 
this, one first has to diagonalize the Bogoliubov Hamil¬ 
tonian. Importantly, one cannot diagonalize the Bogoli¬ 
ubov Hamiltonian by a simple unitary transformation as 
the resulting quasiparticles would not fulfill bosonic com¬ 
mutation relations. In contrast, the diagonalization has 
to be performed in a symplectic manner. Consequently, 
one cannot apply the definitions of topological invariants 
for usual noninteracting systems directly, but has to re¬ 
spect the symplectic nature of the diagonalization proce¬ 
dure. In the next section, we briefly recall this procedure. 
Then we can define a symplectic extension of the macro¬ 
scopic polarization constituting a topological invariant 
that characterizes the Bogoliubov-Bloch bands. 


This leads to the appearance of an effective chemical po¬ 
tential in H^ b \ which reads 


Meff — jyj- (C^mCum Cf,mCt,m) 


^0 (Cf,raCt,ra+l CpmCbm+l 4" C.C.) 

— ^so ( — C^mC-bm+l 4" C|%mC|o 7 ™+1 T C.C.) 


+ 2p^ x ^C im O im C- ,mCs' ,rm 


( 7 ) 


where Ct/i,m denotes now the stationary point of the 
Gross-Pitaevskii function Eqp. 


C. Bogoliubov Hamiltonian in momentum space 

We proceed to work in phase II where (Ct,m,Cbm) = 
(1,0). As we have there a translational-invariant Bo¬ 
goliubov Hamiltonian, we can perform a Fourier trans¬ 
formation and obtain a Bogoliubov Hamiltonian of the 

form m 



III. BOGOLIUBOV THEORY 

A generic Bogoliubov Hamiltonian can be written in 
the form = \ (a* - , a) H (a, a^) T , where 

/ H (“) H(« \ 

h_ V[ h(/3 T (12) 

and at = ^d ^ =1 • • • dt =jV -^. The label p may denote the 

position, momentum, or internal degree of spinor bosons. 
The matrix is Hermitian and is symmetric. We 
also assume that H is positive definite. 

This Hamiltonian can be diagonalized with the ansatz 
(at, a) = ^£>t,b^ Tt, where T denotes a 2A f x 2J\f pa¬ 
raunitary matrix nansnns]. The new operators b shall 
also fulfill bosonic commutation relations. To this end, 
one has to require that 

T^cr^T = ct z Tcr^Tt = cr z , (13) 

with the diagonal matrix ( cr z ) l v = Sij'cri and 07 = 1 for 
l < A f and 07 = — 1 otherwise. After inserting the ansatz 
into the Bogoliubov Hamiltonian, one easily sees, that 
the Hamiltonian is diagonalized if 

HT = a z T^ E _ E ) , (14) 

where E denotes a diagonal matrix E = 
diag [Ei,..., Ej\f] with Ei > 0. As a consequence 
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of Eq. (14) and of the symmetric structure of H, the 


paraunitary matrix T can be written in the form 


T = 


U V* 
V u* 


(15) 


hnv = 0 and the boundary of the Brillouin zone /q nv = tt 
are invariant under inversion as — /q nv mod 2i r = /q nv . 


For these momenta relation (16) exhibit a strict con¬ 
straint for the Hamiltonian H^. . 


with U and V being J\f x J\f matrices. In the following 
we denote by U (V) the particle (hole) part of the ex¬ 
citations. Consequently, only the first M columns of T 
contain independent solutions. The other Af columns re¬ 
semble exactly the same Hamiltonian as the first one, also 
having a positive energy. Thus, there are only positive 
excitation energies. 

Finally we remark that the Hamiltonian in momentum 
space H/c in Eq. <§ does not necessarily have the form 
(12). This problem can be solved by formally combining 
the entries of H& and in an enlarged matrix. 


IV. TOPOLOGY OF BOGOLIUBOV 
EXCITATIONS 

A. Symmetry considerations 

Let us assume that there are symmetries S j transform¬ 
ing a noninteracting Hamiltonian as SjH^Sj 1 = 
with a,/3 = ±1. Depending on the value of a,/?, the 
properties of S j and the dimension of the system one 
finds different topological classes [23 EH- For exam¬ 
ple, for a = —/3 = 1 and S j = V = V~ 1 = V\ the 
single-particle Hamiltonian is invariant under inversion, 
which is the focus of our article. For the noninteracting 
Hamiltonian (|9| the inversion operator can be written as 
V — cr z , where a z denotes the usual 2 x 2 Pauli matrix. 

Next we turn our attention to the symmetry relations 
of H&. Thereby, our approach is to consider the symme¬ 
try operator S B = 12 ®Sj, where S j denotes a symmetry 

of and I 2 is the 2 x 2 identity matrix. 

It is natural that the symmetry of the Bogoliubov 
Hamiltonian has a block structure as otherwise the sym¬ 
metry would relate a pure particlelike excitation with a 
particle-hole-like one. Due to the appearance of the ad¬ 
ditional matrices H^ 1,2 \ the symmetries S j of do 
not necessarily create symmetries of H&, so this sym¬ 
metry can be lost. Consequently, the interaction of the 
particles can change the topological classification. 

However, let us assume that the Bogoliubov Hamilto¬ 
nian and the symplectic extension of the inversion sym¬ 
metry V B = I 2 ® V fulfill 


V B U k V B =H_ fc . (16) 

This condition is fulfilled for the Hamiltonian (j 8 |. For 
most k values, this does not impose a constraint as the 
symmetry just connects the Hamiltonian at k with that 
at —k. Due to the periodicity in position space, the mo¬ 
mentum is only defined within the Brillouin zone which 
we assume to have length 27r. Therefore, the momentum 


B. Symplectic Polarization 

For a single-particle Hamiltonian the so-called macro¬ 
scopic polarization constitutes a topological invariant [23j 
1124] . We now want to formulate a symplectic general¬ 
ization of the macroscopic polarization determining the 
topology. 

For the following derivations we consider systems with 
a discrete basis. The extension to continuous systems 
works analogously, but one has to be careful with the 
dimension of the basis. 

Let T be the solution of Eq. ( p~4| ) in position space 
with periodic boundary conditions. The rows of T can 
be labeled with (m, /), where m E {—M/2 + 1,..., M/2} 
denotes the position of the unit cells and l an internal 
degree of freedom within the unit cell. The columns of 
T contain the eigenstates of ( fl4| ). As we consider peri¬ 
odic systems, they can be labeled with the indices (fc, A), 
where k = 2tt rh/M with rh E {—M/2 + 1,_, M/2} de¬ 

notes a quasimomentum within the Brillouin zone and 
A E {1,..., 2 £} denotes the band index. Let us further 
denote the (fc, A)th column of U (V) by T™ k (T” k A ^) as 
being (A f = M x £)-dimensional vectors. 

Their entries can be expressed in terms of Bloch func- 

tions ( T (U) W) = VM eikmt Ux),v where c e {«,«}• 

The periodic parts are the solutions of the Bogoliubov 
equation in momentum space 


Hfct fc = cr z t k ^ Efe _ E ^ ) , (17) 

where H/e is the matrix in (| 8 |. The paraunitary matrix 
tfc has dimension 2 C x 2 C. It also fulfills t\cr z tk = cr z 
and tkcr z t\ = cr z . More precisely, for A < C the re¬ 
lation reads = ^, A ),z mod/:’ with c = u for 

l < C and c = v otherwise. For A > C we have 
(tfe)z,A = mod /:)*’ with c = v for l < C and 

c = u otherwise m- 

Analogously to the noninteracting case we define the 
corresponding Wannier functions for M —>■ 00 as 

w \,m,l = ^ J bz dkelkmt lk,\),l > ( 18 ) 


where BZ denotes the Brilloin zone. Here we explicitly 
distinguish between particle c = u and hole c = v con¬ 
tributions to the Wannier function. Before defining the 
polarization, we have to sort the bands A. As can be seen 
in Eq. (14), we always have pairs of energies ±E. We con¬ 
sider here only the columns with A < C corresponding to 
positive E. Due to (15), the columns corresponding to 
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negative energies contain only copies of A < C. We con¬ 
sider the bands up to an energy E max which shall be in 
a band gap. Sorting the columns with E > 0 by energy, 
we denote the band with the largest energy E\ < E x n ax 
with A max . Thus we consider the bands A < A max so 
that there is an energy gap between A max and A max + 1 . 
For the Hamiltonian we depict some dispersion rela¬ 
tions in Fig. [ljd) . As C = 2 , we have two bands and the 
spectrum is gapped between A = 1 and A = 2. 

For c = u,v we separately define the corresponding 
contributions to the macroscopic polarization to be 

M/2 

Pc ^m“oo E E ( W \,m,lY (!9) 

TYl— — M/2 + 1 A < A max 

With these definitions we can define the symplectic 
polarization as the difference of the particle and hole po¬ 
larization contributions 


where B k denotes the sewing matrix. Importantly, it 
can only mix states being degenerate. When there are 
no degeneracies, the sewing matrix reduces to a diagonal 
matrix with elements of unit modulus. As we assumed 
that our system is gapped between A max and A max + 1 , 
the sewing matrix B k has a block-diagonal structure. We 
denote the block referring to the band below the gap by 
B <? fc. For the inversion-invariant momenta /q nv , its de¬ 
terminant is a product of the eigenvalues of V B regarding 
the eigenstates, V B t kiavt \ = V\(hnv)t kinv ,\, thus 

det [B <jfcinv ] = P[ r+fcinv). (24) 

A<A m ax 

The sewing matrix obeys the same transformation 
rules as tfc, 

B\a z B k = t\V B cr z t-kCr z t^_ k cr z V B tk 

= a z . (25) 


p = p — p — _ 

1 s — 1 U 1 V - ~ 


T [ dkA(k), 

T7T J bz 


j bz 

where we introduced the Berry potential 


A(k) = i Tr 


A<A n 


+ id 

V X CT z t k CT z ( m t k 


( 20 ) 


( 21 ) 


We define the matrix (r^) ■ = djj'dj,x as being a 2 Cx 

2 C matrix. The symplectic polarization of the bands 
A < A max determines, whether or not there is an edge 
state between the bands A max and A max + 1 . 

For the noninteracting case the symplectic polarization 
reduces to P s -+ P u and coincides with the macroscopic 
polarization of Ref. [24]. Equation ( 21 ) agrees with the 
Berry potential of Ref. [16. found in the context of a 
bosonic Chern number. Yet, in that article there is no in¬ 


terpretation in terms of the symplectic polarization ( 20 ). 
A proof of the last step in (20) is given in Appendix |B 1 


We also prove in Appendix |B 2 | that the Berry potential 
is real valued. 


Using the sewing matrix, we link the symplectic polar¬ 
ization P s to the eigenvalues of the symplectic inversion 
operator. Analogously to the noninteracting case J23], we 
need to relate the Berry potential at k and — fc, but re¬ 
specting the symplectic structure of the eigenstates. We 
find 


A(-k) = -A(k) + idk In [det (B <j fc)] , (26) 

which we prove in Appendix |B 3| Using this we finally 
find 

P * = T- I" dk[A(k) + A(-k)} (27) 

2^ Jo 

= 2 ": J dkd k In [det (B <>fe )] (28) 

= A- {In [det (B<+] - In [det (B <)0 )]} (29) 

Z7T 

= + ln II ^(°)^(+ • (30) 

A < Amax 


C. Topological invariant 

The symplectic polarization of inversion-invariant sys¬ 
tems is strictly quantized to the values P s = m, \ + m 
with m G Z as the one for noninteracting systems m- 
The proof also works essentially as in the noninteracting 
case, yet one has to respect the symplectic structure of 
the eigenstates. We first define the sewing matrix con¬ 
necting the state at k with the one at —k. If we have a 
solution of the Bogoliubov equation 0 in momentum 
space tfc, then V B t k diagonalizes H_/c. Thus, one can 
connect the paraunitary matrices at k and — k as 

B k = Y_ k a z V B t k <^> ( 22 ) 

t_fc = V B t k <r z Bl, (23) 


In the last step we have used that the eigenvalues of 
the inversion operator are r]\(k[ nv ) = =bl. Represent¬ 
ing the eigenvalues in the form r]\(k[ nv ) = 1 = e z2?rm 
or rj x (k inv ) = — 1 = e*( 7r+27rm ) finally proves that P s = 
m, \ + m. 

D. Polarization of the lowest band 

The Bogoliubov excitations of a BEC typically exhibit 
a Goldstone mode in the lowest band denoted here by 
A = 1. This means a linear dispersion for small \k\ and 
thus Ek ,i oc |fc|, which can be seen in Fig. 0d). The 
solution at k = 0 resembles the mean-field solution To 
obtained by the Gross-Pitaevskii equation in the form 
to,i = (^o,i,U),i) T = (To,— To) T . Yet, this solution is 
not normalizable according to as tg z t o,i — 0 . 
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The fact that to 1 is not normalizable is an obstruction 


for defining the Berry potential in Eq. (21) at k = 0. Nev¬ 
ertheless, we argue here how to circumvent this obstacle. 
We use a slightly modified definition for the symplectic 
polarization to respect the difficulties of the lowest band, 
namely, 


P s = — lim 
27r si o 


/ —(5 

dkA(k) + / dkA(k) 

-7T J6 


(31) 


Of course, relation (26) for k ^ 0 is not affected and 
therefore one can adopt the derivation up to line (29) by 
including a limit operation so that 


Ps = f {in [det (B*)] - lim In [det (B*)] ) . (32) 

Z7T <40 J 

The crucial point is to discuss the limit. We also assume 
that the lowest band is non-degenerate in a finite region 
around k = 0. Consequently, the reduced sewing matrix 
has the form 

B < ‘ = ( ^ B<,„ ) ' (33) 


The submatrix B <? / e behaves regularly for k —>> 0 and 
does not cause any problems. So we just have to dis¬ 
cuss the implications of To this end, for a given 

Gross-Pitaevskii solution To, we define a normalization 
function /*. > 0 so that 

lim f k t k , i = (¥o, -^o) T • (34) 

fe->- 0 

The exact shape of fk is not crucial for our discussion, 
but in agreement with inversion symmetry we require 
fk = f-k • For k 0, the relation between the solutions 
at momenta connected by inversion reads, according to 

(HI), 

t-M = e ib ^V B t kA . (35) 


This relation is not well defined at k = 0. Therefore, 
we multiply fk so that the limit k —» 0 exists on both 
sides of the equation. Consequently, 

e ibo,ipB = 1 (36) 


and this constrains 6 q 5 i = 27t m or 6 0 ,i = tt(2 m + 1) de¬ 
pending on the eigenvalue of V B regarding (To, —To) . 
Conclusively one can say that, although the limit k 0 
of the state is not well defined, the phase is up to 2tt m, 
so the final outcome in (30) is not affected. 


E. Application 

Here we return to the Hamiltonian and investigate 
its topology. As the Bogoliubov Hamiltonian and the 
matrix r z = I 2 <£) cr z fulfill 


(37) 
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FIG. 2: (Color online) (a) and (b) Spectra of the Bogoli¬ 
ubov excitations with fixed boundary conditions for M — 30 
sites in phases Ha and lib, respectively. Before performing 
the Bogoliubov diagonalization, we first minimize the Gross- 
Pitaevskii functional. The resulting wave function of the con¬ 
densates close to the boundaries is depicted in the insets. Due 
to the fixed boundary condition, the density of the condensate 
is lower at the boundaries. Using the condensate wave func¬ 
tion, we then perform the Bogoliubov diagonalization. The 
excitation energies are all positive and real valued, as we have 
prepared the wave function in its ground state. One clearly 
identifies the midgap states, (c) Wave function of one of these 
midgap states. 


the operator r z is an inversion symmetry. In Fig. [ljd) 
we depict some dispersion relations of the system. 


To determine the topology of our model we have to con¬ 
sider the inversion-invariant momenta /q nv = 0, n. Let us 
consider the lower band. As the mean-field wave func¬ 
tion (Cm,t?Cm,4) — (1,0) is constant for the parameters 
in phase II, the eigenvalue under inversion is rjk=o, 1 = 1. 
The gap closes just at the boundary of the Brillouin zone 
so that the symplectic polarization P s can only change 
there. To investigate this, we consider the eigenstates at 


'T z^A-k^ z — H — fc, 
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k = 7r, which read 

^7r,i = 4 a/ ^o(pX + v$), ^,2 = 2(cj — px + PXtt)> 


^ 77,1 


PX + 2^o - 2^/t' 0 (px + t'o) \ 
0 




PX 

0 


J 


(38) 


5 ^ 77,2 - 


(39) 


Obviously, they fulfill T z t^^\ = — (— l) x t 7T ,\. Conse¬ 
quently, the topological invariant (30) changes at the de¬ 
generacy point E , 7r? i = Thus, the boundary be¬ 

tween the topological phases is 


^o,tpt — 

—PX + y / 2(px) 2 - 2px(^ + PXn) + (^ + PXtl) 2 ^ 

We depict this topological phase boundary also in the 
phase diagram in Fig. [ljb) for xtt = 0- One can see > 
that the product p\ has a strong impact on the topology 
of the system. The topological invariant for the lower 
band is P s = 0 for z/o < ^o,tpt and changes to P s = \ for 
> z^o,tpt • Accordingly, the system is in a topologically 
trivial phase or a nontrivial phase, respectively. 


F. Edge states 


Although we have defined a topological invariant, there 
is still the question about the physical consequences of it. 
In contrast to fermionic systems, where the polarization 
is an actual physical quantity, in bosonic systems this is 
not the case as not all momenta of a band are equally 
occupied. However, the symplectic polarization (20) de¬ 


termines the existence of edge modes in a system with 
finite length and fixed-boundary conditions. 

As a demonstration, we consider our instructive model 
in Fig. [2] For this illustration, we do not assume an ad¬ 
ditional harmonic confining potential. Due to its topo¬ 
logical origin, the edge states are robust in the presence 
of moderate perturbations 39 . We also refer to Ref. [40 
for the creation of sharp boundaries. 

Importantly, first we have to determine the mean-field 
wave function by minimizing the Gross-Pitaevskii func¬ 
tional. Here the Gross-Pitaevskii mean field has no uni¬ 
form density at all sites due to the boundaries. We depict 
the part of the condensate close to the left boundary in 
the insets of Figs. [2^a) and (b). One can see that the 
density is smaller close to the edges. Away from the 
boundaries, the mean field is approximately constant, so 
the results derived in Sec. m are still valid. 

On top of the mean-field wave function, we perform 
a Bogoliubov diagonalization in position space and de¬ 
pict its spectrum in Figs. [2ja) and (b). In Fig. |2|a) the 


system is in the trivial phase Ha, so no midgap states 
appear. In contrast, one clearly identifies two midgap 
states within the spectrum in Fig. §b) depicting the 
spectrum for parameters in phase lib. In Fig. §c) we 
plot the wave function of one of these states. We find 
that it is strongly localized on the edges. The interpreta¬ 
tion of these edge states works analogously to that in the 
fermionic case ESI. There, each edge state contributes 
half an electron to each boundary. Thus, one particle 
splits into two half particles. In the bosonic case, corre¬ 
spondingly, each edge mode can be considered to consist 
of two half modes at the boundaries. 


V. CONCLUSIONS 


We investigated the topology of Bogoliubov excitations 
in inversion-invariant systems of interacting bosons. To 
characterize the topology, we extended the definition of 
the macroscopic polarization in a symplectic manner. We 
called this quantity symplectic polarization which is de¬ 
fined in Eq. (20) as the difference of the particlelike and 
holelike polarization contributions. As in noninteracting 
systems with inversion symmetry, this quantity can be 
expressed by the inversion eigenvalues of the states at 
inversion-invariant momenta. In an instructive example 
we showed, that the topological invariant strongly de¬ 
pends on the condensate density, so the interaction be¬ 
tween the particles modifies the topology of the excita¬ 
tions. 


The definition of the symplectic polarization can also 
be applied to analyze inversion-invariant systems in 
higher dimensions. In this case we expect that an invari¬ 
ant defined as the product of the inversion eigenvalues 
of the states at the inversion-invariant momenta predicts 
edge states [23]. Furthermore, one can link also the sym¬ 
plectic polarization of the one-dimensional system to the 
Chern number in two dimensions Another possible 
application is in inversion-invariant systems consisting of 
arrays of (pseudo)spins with large angular momentum 
such as in Refs. mmm, where fluctuations above 
the mean-field treatment are diagonalized by Bogoliubov 
theory. 

A drawback of the bosonic edge modes in the excita¬ 
tion spectrum is that they are weakly occupied. How¬ 
ever, this obstacle could be circumvented by engineering 
a bosonic system in a driven setup within Floquet theory 
analogously to [SSI . 

Importantly, the symplectic polarization discussed 
here can be used to define the symplectic generalization 
of the time-reversal polarization [44] . This can be used to 
analyze the topology of time-reversal invariant systems of 
interacting bosons. 
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Appendix A: Minimization of the Gross-Pitaevskii 
functional 

The Gross-Pitaevskii functional reads 


M/2 


F’gP = ^ Eo,m + FV,m, 

(Al) 

m = —M/2+1 



where 


U (CpmCUm — C+mCt,m) 


(C+mCt,m+l — CbmCUm+l + C * C *) 


Uso (C+mC4,ra—1 — C+mC-Rra+l + C.C.) , 

(A2) 

^ ^ PoXs,s'Cs,mCs',mCs,mCs',m- 

(A3) 

s,s'= t,4 
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Motivated by Ref. f33], we use the modified ansatz 

( A’ m ) m. Ci ( COs9 a ) e ikm + C 2 ( COSe Q ) e~ ikm . 
v Ci,m ) \ -1 sin 9 J \ i sin 9 ) 

(A4) 

The variational parameters are Ci, C2, 0, and k. As 
we work at a fixed particle number, we have to respect 
the constraint |Ci| 2 + |C 2 1 2 = 1. Due to this ansatz the 
noninteracting part of the energy functional reads 

m 

= M cos 26 (—w — 2z/ 0 cos k) + M2z/ so sin 26 sin k. (A5) 

Importantly, the noninteracting part does not depend on 
the coefficients C\ and C 2 which reflects the inversion 
symmetry of the system. Accordingly, the interaction 
terms turn out to be 

F Ev,m = Mpo (1 + 2/3) (xtt cos 4 9 + xw sin4 A 

rri 

+ Mp 0 2xu (1 — 2/3) cos 2 6 sin 2 6 , (A6) 

where we define (3 = \C \\ 2 |d7 2 1 2 with /3 G (0, |). We im¬ 
mediately see that Eqp is a linear function of /3. There¬ 
fore, the minimum can only be located at /3 = 0 or 

P = 1/4- 

Next we derive a relation between k and 6. To this end 
we take the derivative of Eqp with respect to k. After a 
short algebraic manipulation, we obtain 

tan 20 =-— tanfc. (A7) 

"so 

Having done these analytical preparations, we can 
numerically minimize the Gross-Pitaevskii functional, 
which is now just a function of essentially one variable, 


thus Eqp = Eqp(6,(3) as /3 G { 0 , The result is de¬ 
picted in Fig.[ljb). For comparison, we also directly min¬ 
imize the Gross-Pitaevskii functional numerically, where 
we also found the localized ground-state wave function 
for xPo < 0. 


Appendix B: Berry potential 


1. Details of the derivation 


We start with the final expression of the Berry poten¬ 
tial and perform the proof from the end. First, we insert 
the representation (15) for tk and perform the multipli¬ 
cations with cr z so that we obtain 


A(k) = i F Tr 

T\(T z tl(T z (<9fctfc) 

A<A max 


= i E Tr [ r d 

A<A max 

< 4 -ytu 

V -vl A; n T _ k ) dk 

= * E (4+ u fc - v l+ v fc) A A - 


A<A max 


(Bl) 


u fc v -/c V 
v * u -/c A 

(B2) 

(B3) 


We evaluate the matrix product by inserting a complete 
1 of the basis states l of the unit cell so that 

m=iJ2 E ( B4 ) 

l A < Vax 
C=U,V 


using a(u ) = +1 and a(v) = —1. We have also used that 
the columns of and v/~ are the periodic parts of the 
Bloch function, namely, t™ k x ^ L and t” k x ^ v respectively. 
We continue by inserting unity so that 


lL dkm =iL dkdk '^ + ( 

c=u,v 

= lim [ dkdk' 

(2tt) 2 M^oo7 bz A a< A 

t A< A lrl ax TYl —- 
C=U,1) 

= — \ lim [ dkdk'^S^ 'S 

(2 % £ 


M/2 

E < 

■M/2+1 


= ~ tv xo hm 
(27r) 2 M-)>oo t 

— P — P 

1 U 1 V 


' '■ A max 

C=U,V 


^ ''max 
C=U,V 

I 


)+ +,A),/ 

(B5) 


(B6) 


(B7) 


(B8) 


(B9) 


This finally proves Eq. (20). 
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2. Real valueness of the Berry potential 


To prove that the Berry potential (21) is real valued, 
we calculate 

A*(k) = -i ]T {Tr \r x a z t\cT z (3 fc t fe )] }* (BIO) 

= -i ^2 Tr VztkCTz r 

A^ A max 

= -i ^2 Tr (^fetfe) ^ z t k cr z T 

A^Amax 

= i E Tr t£<r z (d k t k ) cr z T 

A (0 A max 

= A(fc). 


(BH) 

(B12) 

(B13) 

(B14) 


3. Symmetry relation 


Here we prove relation (26). By definition we have 


A(-k) = i E Tr r Acr z t^ fc cr z (<9_ fc t_ k ) (B15) 

A ^ A ma x 

= *^Tr T x B k <j z t\V B <J z (d-kP B tk<J z B^ (B16) 

A 

= i J2 Tr [TxB k a z tlV B a z V B (c>_ fe t fe ) cr z B*] (B17) 

A 

+ i ^ IV \r\B k <T z t ] k V B (T z V B t k <T z (B18) 

A 

=-A(k) - i ^2 Tr T x a z B k a z d k Bl . (B19) 

A 

= -A(k) - iTr [B <)fc 4B t <)fe ] (B20) 

Next we recognize that 


cr z B k cr z Bl = 1 


(B21) 


due to Eq. (|25|). As the sewing matrix has a block diag- 

(B22) 


onal structure, we find 


— 1 . 


This relation means that B <? / c is unitary. We expand it 
in terms of its eigenstates | i) such that 

B<,fe = E e * 10 o ^ B E = E E 10 (0 • ( B23 ) 

i i 


Both the eigenstates and the eigenvalues depend on k. 
For a notational reason we suppress the index in the fol¬ 
lowing. Now we can prove relation (26): 


Tr 




= Tr 


= Tr 


Tr 


Tr 


E e * 10 01 I ^E e i 1 io oi 

i \ j 

E>io oiEOaE) looi 

i 3 

E e *o E 2 e 7 x io>) oi 

i 3 

E ei io (^lET io (<aod 


(B24) 


(B25) 


Evaluating the traces in the eigenbasis of B <? / c we find 


Tr 


B< )fe <9 fc B< ife 


= E^lne* 1 + E^* \ dk 10 + ( d k (*|) \i) 

i i 

= -d k In 11 + ^2 (9k (i\ i)) 

i i 

= -d k In [det (B <)fc )]. (B26) 


In the derivation we have also used that | e* | = 1, so that 
dfcep = ~d k ei. 



































